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Abstract 

The aim of this paper is to dynamicaUy constrain the quadrupole 
mass moment J2 of the main-sequence HD 209458 star. The adopted 
method is the confrontation between the measured orbital period of its 
transiting planet Osiris and a model of it. Osiris is assumed to move 
along an equatorial and circular orbit. Our estimate, for given values of 
the stellar mass M and radius R and by assuming the validity of general 
relativity, is J2 = (3.5 ± 385.1) x 10~^. Previous fiducial evaluations 
based on indirect, spectroscopic measurements yielded J2 ~ 10~^: such 
a value is compatible with our result. 
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1 Introduction 

The quadrupole mass moment J2 is an important astrophysical stellar pa- 
rameter related to the inner structure and dynamics of a star (Paterno et 
al. 1996; Pijpers 1998). 

In this paper we dynamically constrain the quadrupole of the main- 
sequence star HD 209458 from the measured orbital period P of its transiting 
planet HD 209458b 'Osiris' (Charbonneau et al., 2000; Henry et al., 2000). 
It is the best-studied transiting planet to date, mainly due to its proximity 
and the resultant high apparent brightness of its parent star. 

Until now, only a fiducial value J2 = 2 x 10"^ by Miralda-Escude (2002), 
based on spectroscopic measurements of rotational velocity (Queloz et al. 
2000), exists for HD 209458. Also Winn et al. (2005) assumed J2 ~ 10"*^. 
More generally, Miralda-Escude (2002) investigated the possibility of dy- 
namically measuring the quadrupole of a star from the node and periastron 
precessions of a transiting planet which cause a time variation of the dura- 
tion of a transit; the periastron precession also induces a variation of the 



transit period. Such orbital perturbations should be measured from an ac- 
curate photometric analysis of the light-curve of the transiting planet, but, 
until now, such a proposed strategy has not yet been implemented. Winn 
et al. (2005) pointed out that, for J2 ~ 10^^, the quadrupole node preces- 
sion would amount to about 4 arcseconds per year; measuring such an effect 
would require high-precision photometry spanning several years (Winn et 
al. 2005). 



2 The use of the orbital period of Osiris 

The Newtonian gravitational potential U of an oblate star of mass M and 
equatorial radius R can be written as 

^ = + ^^(3cos^^ - 1), (1) 

where 9 is the co-latitude angle. For the orbital period of a planet of mass 
m in a circular and equatorial {9 = 7r/2) orbit of radius a eq. (1) yields 



P(N) ^ p(0) + p(J.) ^ 2J—-^^ - ^ ^^^'-^^ ■ (2) 

In fact, in addition to eq. (2), there is also a post- Newtonian, general 
relativistic part (Soffcl 1989; Mashhoon et al. 2001; lorio 2005; 2006) to be 
added; for circular orbits and m <C M (see Section A) it is 



p(PN) _ ^^r^/G{M + m)a 



(3) 



where c is the speed of light in vacuum, so that 

p ^ p(0) _^ p(J2) + p(PN)_ (^4) 

In the case of Osiris, eq. (2)-eq. (3) yield a reliable model of its orbital period 
because the eccentricity e was recently evaluated to be e = 0.014 it 0.009 
(Laughlin et al. 2005) and the inclination angle of the orbital plane to 
the star's equator, determined by means of the Rossiter-McLaughlin effect 
(Rossiter 1924; McLaughlin 1924), should not be larger than about 5 deg, 
(Winn et al. 2005). Moreover, there are currently no observational evidences 
of the presence of other bodies around HD 209458 (Brown et al. 2001; CroU 
et al. 2005; Laughlin et al. 2005; Agol & Steffen 2006) which may require 
the introduction of additional perturbing terms in eq. (4); the inclusion of 
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the general relativistic correction of eq. (3) is required because it amounts to 

about 0.1 s, while the errors in the most recent measurements of the Osiris' 
orbital period are 0.016 s (Wittenmyer et al. 2005) and 0.033 s (Knutson et 
al. 2006). 

Thus, the HD 209458 quadrupole mass moment can be determined by 
comparing the model of eq. (2)-eq. (4) to the measured period p{«icas)^ 
determined in a purely phenomenologically way from combined photomet- 
ric transit and spectroscopic radial velocity techniques, independent of any 
gravitation theory, and solving for J2 

2p{meas) ^ ^ 4/a\2 2G{M + m)a 

By using eq. (5) and the system parameters derived for M = l.OTM© and 
R = 1.137R0 and p("^e^) = 3.52474554 d (Wittenmyer et al. 2005), we 
obtain 

J2 = 3.5 X 10-^ (6) 

It must be noted that the obtained result is free from any a priori, 
'imprinting' effect by J2 itself. Indeed, M and R are kept fixed, and a, 
determined from^ 

Kf^P^ o?rrr' sin'^ i 

" (m + M)3 ^ ' 

which is independent of any model of the orbital period, is not affected by 
J2 over timescales longer than one full orbital revolution. 
Let us now evaluate the uncertainty in J2 as 

(5 J2 < <5 jf ) + + Jji'"^ (8) 

For the same values of M and R as before and (5p("ieas) ^ q gig g (Witten- 
myer et al. 2005) we have 



5j!p^ = 



8 a I pi^ee-s) l G{M+m) . 2G{M+m) 



P(°'«^°) / Ga I 2Ga 



Sa = 2.404 X 10"^, 
6m = 1.437 X 10-^, (9) 



2^yG{M+m)a 
3^ 



^p(meas) ^ 5 ^ iq~6 _ 



^Km = (2l£M.^ sin j — ^_^L_ j irw') sini is the projected semiamplitude of the star's 
radial velocity. 
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The stellar mass was not included in the least-square solution by Wittenmyer 
ct al. (2005) also because its determination is more model-dependent than 
the other parameters. The range of allowable values 1.06 it 0.13 solar masses 
(Cody &: Sasselov 2002) was, instead, used; it comes from observational 
errors in temperature, luminosity, and metallicity as well as systematic errors 
in convection mixing-length and helium abundance. The resulting scattering 
in the determined values of J2 is 

3.2 X IQ-^ < J2 < 3.7 X IQ-^ (10) 

Thus, we can state that the model-dependence of our estimate amounts to 
5 X 10~^, so that a conservative estimate of the total uncertainty in J2 is 

5J2 < 3.851 X 10"^. (11) 

3 Discussion and conclusions 

In this paper wc dynamically constrained the quadrupole mass moment J2 
of the HD 209458 star from the orbital period P of its transiting planet 
Osiris, assumed to be in a circular and equatorial orbit. Its measured 
value— determined in a phenomenological way, independent of any gravi- 
tational theory— was compared to an analytical model including the New- 
tonian part, constituted by the usual Keplerian component and the term 
induced by J2, and the post-Newtonian, general relativistic correction. The 
inclusion of the latter term, which is of the order of 0.1 s, is motivated by 
the ~ 0.01 s level of accuracy reached in measuring P. By keeping the stel- 
lar mass M and radius R fixed to values within a range determined from 
stellar evolution models and temperature/luminosity measurements, by as- 
suming that general relativity is valid in the HD 209458 system as well and 
that Osiris is the only planet affecting the motion of its parent star in a 
detectable way, we obtain J2 = (3.5 it 385.1) x 10~^. While the uncertainty 
due to the error in the orbital period amounts to ~ 10~^ only, the Osiris' 
mass and semimajor axis boost the total bias to ~ 10~^. Previous fiducial 
evaluations based on indirect, spectroscopic measurements giving J2 ~ 10~^ 
are compatible with our result. 

In order to make easier a comparison with our results, in Table 1 we quote 
the numerical values used for the relevant constants entering the calculation. 
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Table 1: Values used for the defining, primary and derived constants 
(http:/ / ssd.jpl.nasa.gov/?constants#ref ) . 



constant 


numerical value 


units 


reference 


c 


299792458 


m 


(Mohr & Taylor 2005) 


GMq 


1.32712440018 x 10^° 




(Standish 1995) 


G 


(6.6742 ± 0.0010) x lO'^^ 


kg^i m'^ s^^ 


(Mohr & Taylor 2005) 




6.95508 X 10« 


m 


(Brown & C.-Dalsgaard 1998) 


1 mean sidereal day 


86164.09054 


s 


(Standish 1995) 



A The post-Newtonian correction to the orbital 
period 

In fact, the post-Newtonian gravito-electric correction to the orbital period 
does depend on both the eccentricity e and the initial value of the true 
anomaly /o according to 



p(PN) 



Svr 



^/G{M + m)a 



F{e,fo), 



with (Soffel 1989; Mashhoon et al. 2001) 



FieJo) = 3-- 



and 



3 (l + ecos/o)2' 
mM 



(12) 



(13) 



(14) 



(M + m)2 ■ 

Depending on e and v, F-and P(^^) -vanishes for those values of /o which 
satisfy the relation 



cos /o = - 

e 



- 1 



(15) 



it may also happen that the absolute value of the right-hand-side of eq. (15) 
is larger than 1, so that P^^^) ^ 0. 

In the case of the HD 209458 system, v ~ 10~^; recent refinements in 
the Osiris ephemerides yields e < 0.023 (Laughlin et al. 2005), so that 



9-1/ 



(16) 
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F never vanishes, ranging from 0.90 to 1.09. Thus, the difference between 
the maximum and the minimum values of p(P^) is 0.019 s at the most: it 
just lies at the edge of the precision with which the orbital period is known, 
i.e. 0.016 s (Wittenmyer et al. 2005) and 0.033 s (Knutson et al. 2006), so 
that we can approximate F to unity. 
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